Origin of the large thermoelectric power in oxygen-variable i?BaCo205+a; {R=Gd, Nd) 
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Thermoelectric properties of GdBaCo2 05+a; and NdBaCo2 05+a; single crystals have been studied 
upon continuous doping of C0O2 planes with either electrons or holes. The thermoelectric response 
and the resistivity behavior reveal a hopping character of the transport in both compounds, providing 
the basis for understanding the recently found remarkable divergence of the Seebeck coefficient at 
X = 0.5. The doping dependence of the thermoelectric power evinces that the configurational 
entropy of charge carriers, enhanced by their spin and orbital degeneracy, plays a key role in the 
origin of the large thermoelectric response in these correlated oxides. 

PACS numbers: 72.20.Pa, 72.80.Ga, 72.20.Ee 



Materials with strong electron-electron interactions, 
called strongly correlated electron systems, have recently 
attracted a great deal of attention as a promising al- 
ternative to conventional semiconductors in the field of 
thermoelectric power generation. The hope to find com- 
pounds with superior thermoelectric properties among 
correlated materials is based on the idea that in cor- 
related systems, where the spin and orbital degrees of 
freedom play an important role in transport properties, 
the thermoelectric response can be enhanced by a large 
spin-orbital degeneracy of charge carriers P, ■ 

Naa;Co02 and several misfit-layered cobaltites with 
triangular-lattice C0O2 planes — materials demonstrat- 
ing high thermoelectric response and metallic behavior 
1^ — have been considered as possible examples of com- 
pounds where the spin-orbital degeneracy plays a dom- 
inant role in enhancing the thermoelectric power 0, Q ■ 
However, a conventional Boltzmann-transport approach 
can also explain the coexistence of a large Seebeck coef- 
ficient and metallic conductivity in these compounds , 
blurring the role and relevance of strong electron corre- 
lations. 

Recently, layered cobaltites i?BaCo205+2: (R is a rare- 
earth element) with square-lattice C0O2 planes came into 
focus because of their remarkable transport and magnetic 
properties 0, 0, • ^ variety of spin and orbital states 
are available in i?BaCo205+2; owing to its very rich phase 
diagram 0, where x — 0.5 is the parent compound with 
all cobalt ions in the 3+ valence state. The variability of 
oxygen content in i?BaCo205+a; allows one to dope con- 
tinuously the C0O2 planes with either electrons (Co^"*" 
states) or holes (Co^+ states) and to measure a precise 
doping dependence of the Seebeck coefficient in one and 
the same crystal. Apart from this filling control, one 
can also employ a bandwidth control of transport proper- 
ties using rare-earth elements R with different ionic radii 
[Fig. 1(a)] to clarify the role of band structure parame- 
ters in determining the thermoelectric power in these ox- 
ides. Thus, i?BaCo205-|-2; can provide a suitable ground 
for elucidating the relation between the spin-orbital de- 
generacy and thermoelectric properties in correlated ma- 



terials. 

Here we present a comparative study of the Seebeck 
coefficient and the resistivity in GdBaCo2054.a; (GBCO) 
and NdBaCo205+2; (NBCO) single crystals over a wide 
range of electron and hole doping. Despite the difference 
in the lattice parameters, both compounds show virtu- 
ally the same behavior of the Seebeck coefficient, imply- 
ing that the thermoelectric power in these compounds is 
governed by correlated hopping of charge carriers rather 
than the band structure parameters. An analysis of the 
doping dependence of the Seebeck coefficient, showing a 
spectacular divergence at cc = 0.5, provides evidence that 
the entropy contribution of charge carriers is the main 
source of the large thermoelectric power in i?BaCo205+a;. 

High-quality GdBaCo205+a: and NdBaCo205+2; single 
crystals were grown using the floating-zone technique and 
their oxygen content was modified by high-temperature 
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FIG. 1: (Color online) (a) X-ray Bragg peaks (200), (040), and 
(004), measured in GdBaCo2 05.5 and NdBaCo2 05.5 crystals 
at room temperature (each peak has Cu Ka\ and Cu Ka2 
contributions to the diffraction pattern). Note that the unit 
cells are doubled along the b and c axes. A sketch of the 
crystal structure of iiBaCo2 05.5 is shown in the upper part 
of the panel, (b) The evolution of the room-temperature lat- 
tice parameters in GdBaCo2 05+a; and NdBaCo2 05+a; under 
variation of the oxygen content x. 
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FIG. 2: (Color online) Temperature dependences of the 
in-plane resistivity p{T) of (a) GdBaCo2 05+a; and (b) 
NdBaCo205+a; crystals. 
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FIG. 3: (Color online) Temperature dependences of the 
Seebeck coefficient QiT) of (a) GdBaCo2 05+i, and (b) 
NdBaCo2 05+a; crystals. 



annealing treatments with subsequent quenching to room 
temperature 's']. Resistivity measurements were carried 
out by a standard ac four-probe method. The thermo- 
electric power was measured in a slowly oscillating ther- 
mal gradient of ~ 1 K along the ah plane. The contri- 
bution from the gold wires (~ 2 yN jY^) used as output 
leads was subtracted. 

In parent i?BaCo205.5 compounds, the crystal lattice 
consists of equal numbers of CoOg octahedra and C0O5 
square pyramids [schematically shown in the inset of Fig. 
1(a)]. Upon changing the oxygen content, some oxygen 
ions are inserted into or removed from the RQx planes, 
which changes the number of CoOg octahedra and C0O5 
pyramids and also creates electrons or holes in C0O2 
planes. Figure 1(b) shows the room-temperature lattice 
parameters in the entire available range of oxygen con- 
centrations. Regardless of the oxygen content, the unit 
cell parameters in NBCO are found to be noticably larger 
than in GBCO, in agreement with the larger ionic radius 
of Nd. This difference in the unit-cell size implies the 
difference in Co— O distances and/or O— Co— O angles, 
which determine the one-electron bandwidth; therefore, 
the substitution of the rare-earth element should affect 
the thermoelectric properties, if they are governed by the 
band structure parameters. 

Figure 2 shows the temperature dependences of the 
resistivity in GBCO and NBCO for several oxygen con- 
centrations. Both compounds demonstrate similar piT) 
behavior, which strongly depends on the oxygen content. 
For oxygen concentrations close to the parent composi- 
tion X — 0.5, p(T) curves show a sharp metal-insulator 
(MI) transition upon cooling below approximately the 
same temperature Tmi ~ 360 K, while the MI transition 
seems to be smeared for x away from 0.5. Nevertheless, 
the low-temperature resistivity exhibits a hopping char- 
acter for the entire range < a; < 0.7, which allows us 



to consider the charge transport in terms of hopping of 
localized electrons (Co^+) or holes (Co''+). In the insu- 
lating regime, the resistivity in both compounds quickly 
decreases with hole doping, but remains unchanged or 
even increases with electron doping. This doping asym- 
metry can be explained by the different hopping prob- 
ability of localized electrons (Co^+) and holes (Co^+) 
moving in the background of Co'^"'" ions because of the 
spin blockade of the electron transport 0. 

A comparison of the temperature dependences of the 
Seebeck coefficient Q(T) in GBCO and NBCO for sev- 
eral oxygen concentrations reveals a remarkable simi- 
larity of their thermoelectric properties (Fig. 3). Not 
only the temperature dependences are similar, but also 
the absolute values of Q(T) are virtually the same in 
GBCO and NBCO. At high temperatures, both com- 
pounds show a small, negative, and almost temperature- 
independent and doping-independent Seebeck coefficient, 
which is quite natural for a metallic state Q . On the in- 
sulating side, on the other hand, the temperature depen- 
dences of the Seebeck coefficient are rather complicated 
and do not follow a simple ~ 1 jT law expected for insu- 
lators. (3(r) depends strongly on the oxygen content, be- 
ing negative for electrons and positive for holes, as shown 
in Fig. 4(a) for T = 100 K. Note that the absolute value 
of the Seebeck coefficient decreases rapidly upon doping 
with electrons or holes, in contrast to the asymmetric 
doping dependence of the conductivity, shown in Fig. 
4(b) for the same temperature. The most striking fea- 
ture in the doping dependence of the Seebeck coefficient 
Q{x) [Fig. 4(a)] is its remarkable divergence x — 0.5, 
where it reaches a large absolute value. 

In order to clarify the origin of this large thermoelectric 
response, it is useful to consider what one would expect 
for the doping behavior if i?BaCo205+a; were an ordi- 
nary band-gap material. Figure 4(c) shows the plot of 
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Q versus the logarithm of the conductivity measured at 
T = 100 K - the so cahed Jonker plot 0|. For band 
gap materials, this plot is a curve with a universal shape 
[the pear-shaped curve, shown in Fig. 4(c)], which can be 
scaled by intrinsic parameters related to the energy band 
gap and scattering mechanism. In the extrinsic conduc- 
tion region, the curve Q vs. ln((T) has a slope iks/e, 
where ks is the Boltzmann constant and e is the ab- 
solute value of the electron charge. Upon approaching 
the intrinsic-conduction region with decreasing concen- 
tration of electrons or holes, this slope can only decrease. 
As can be seen in Fig. 4(c), the actual slope of the Jonker 
plot for i?BaCo205+2; increases near the intrinsic conduc- 
tivity region, exceeding the ideal fc^/e slope by several 
times. Moreover, the decrease of Q with electron dop- 
ing does not correlate with the doping dependence of the 
conductivity, which shows no increase because of the spin 
blockade phenomenon 0. Thus, both the temperature 
and doping dependences of the Seebeck coefficient clearly 
deviate from the behavior of conventional band gap ma- 
terials. 

Given the hopping character of the charge transport 
in i?BaCo205+a:, one can reasonably assume that a large 
part of the thermoelectric power comes from the en- 
tropy, which is carried by each electron or hole along 
with a charge. To check this assumption, let us exam- 
ine the doping dependence of the Seebeck coefficient in 
GBCO and NBCO crystals, Q{x), measured at T = 100 
K [Fig. 4(a)]. Although the chosen temperature is pre- 
sumably too low to achieve the true high-temperature 
limit, where the entropy contribution becomes the only 
contribution to the thermoelectric power 0, it should 
be high enough to make reasonable estimations. Note 
that for _RBaCo205+a;, which undergoes a transition to 
a metallic state, the small and doping-independent See- 
beck coefficient at high temperatures is governed by a 
different mechanism jg. 

The entropy contribution AQ to the thermoelectric 
power in i?BaCo205+a; should be proportional to the 
change in the entropy S of the electron system upon in- 
troducing A'^ charge carries at constant internal energy 
E and volume V of the crystal 0, where the entropy 
is determined by the total number of configurations Q, 
i.e., all possible arrangements of introduced carriers in 
the crystal lattice. Hence, AQ is expressed in terms of 
as 
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FIG. 4: (Color online) (a) The doping dependence of 
the Seebeck coefficient Q{x) of GdBaCo2 05+a; (circles) and 
NdBaCo2 05+a; (squares) crystals at T = 100 K. The hatched 
area represents the entire available range of the entropy con- 
tribution to the thermopower AQ„{x) and AQp{x). (b) The 
doping dependence of the conductivity aix) of GdBaCo205+a: 
(circles) and NdBaCo2 05+i (squares) crystals at T = 100 K. 
(c) Jonker plot for i?BaCo2 05+i (see text). 
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Both electron AQ„ and hole AQp contributions are 
shown in Fig. 4(a) by the dashed lines. A comparison 
with experimental data shows that this simple model, 
albeit it gives a rough idea of why the Seebeck coeffi- 
cient can diverge at x = 0.5, is not sufficient for a quan- 
titative description of the doping dependence Q{x) in 
i?BaCo205+2,, especially in the case of electron doping. 

As has been recently emphasized by Koshibae and 
Maekawa 0, the spin-orbital degeneracy of charge car- 
riers plays an important role in systems with correlated 
hopping transport. In i?BaCo205+i^, there are two kinds 
of lattice sites for cobalt ions that differ in the local oxy- 
gen environment [see the inset of Fig. 1(a)], namely, x 
octahedra and (1 — x) pyramids, which provide different 
spin-orbital states for introduced electrons (Co^+), holes 
(Co'*"'"), or host Co'^"'" ions. For electron doping, if one 
knows the distribution of introduced electrons (Co^"*") 
between octahedral and pyramidal positions, the total 
number of configurations is determined by 
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kB dlnfl 
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(1) 



Thus, the description of the doping dependence of the 
Seebeck coefficient Q{x) in our case can be reduced to a 
combinatorial problem. 

In the simplest case of nondegenerate carriers, the 
above approach gives the celebrated Heikes formula ll| 



O - N2p N3a N,, 

" — 920 y2p i/3o i/3p 



A^2o! N2p\ Nsol N3p\ 



(3) 



where Nl is the total number of cobalt sites; A'20 and 
(A30 and Nsp) are the numbers of Co^"*" (Co'^"'") ions 
in octahedral and pyramidal positions, respectively; 320, 
<?2p, 53o, <?3p are the degeneracies of these states. An 
analogous expression can be written for hole doping as 
well. 
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At any given doping level in the range of < x < ^ , the 
equilibrium concentrations of Co^^ and Co'^"'' ions in dif- 
ferent positions are determined by the only parameter - a 
probability p to find a Co^"*" ion in an octahedral position. 
For X close to x = 0.5, where the number of doped carries 
is not large, this probability is determined by the law of 
mass action (A^zp • N3o)/iN2o ■ Nsp) = K(T) [ij, which 
couples the equihbrium concentrations of Co^+ and Co^+ 
ions that are competing for octahedral and pyramidal po- 
sitions. This equation can be analytically solved for p, 
yielding 



K+l 
2 



2{K-l)-a-x) 



(4) 



where K{T) is the equilibrium constant, which depends 
only on the temperature and intrinsic characteristics of 
the system, such as the difference in energies of Co^+ 
(and Co'^+) ions in octahedral and pyramidal positions. 

Using Eqs. (1), (3), and (4), the analytical solution 
for the entropy contribution of doped electrons to the 
thermoelectric power can be written as follows: 



AQ„ — — In 

e 



92o 



A,^+'" A^,~" 



(5) 



where m = p+{x — \ )-|f - Upon electron doping, m can 
change its value in the range from -1 to -1-1, depending on 
the intrinsic characteristics of the system. The entropy 
contribution of doped holes to the thermoelectric power 
in the range \ < x <1 can be obtained in the same way. 

As follows from Eq. (5), there are only two factors 
that affect the doping dependence of the Seebeck coeffi- 
cient: The equilibrium constant AT, which determines the 
occupation of octahedral (pyramidal) positions by doped 
electrons, and degeneracies of states gi, which turn out 
to be the main source of the thermoelectric-power en- 
hancement in i?BaCo205+a; near x — 0.5. Note that the 
entropy contribution to the thermoelectric power is ex- 
pected to be insensitive to the unit-cell size, which gives 
a natural account for the similarity of the Seebeck coef- 
ficient Q{T) in GBCO and NBCO. 

It is well known that owing to a comparable strength 
of the exchange interaction and the crystal field splitting, 
Co'^"'" ions can adopt different spin states, depen ding on 
the oxygen environment or even on temperature 'l2j . It 
is believed that in the parent i?BaCo205.5 at low tem- 
perature, the Co^"*" ions adopt the low-spin (LS) state in 
octahedra and the intermediate-spin (IS) state in pyra- 
mids 0, I3- For Co^+ the crystal field is weaker than 
for Go^^ and thus Co^+ adopts the high-spin (HS) state, 
while for Co^+ the crystal field is stronger, favoring the 
LS state. Any state with nonzero spin S has the degen- 
eracy 25* -|- 1. In addition to this spin degeneracy, there 
can be orbital degeneracy as well. For both IIS-Co^+ 



and LS-Co^+ the orbital state is threefold degenerate be- 
cause of a hole residing in the t2g orbitals. In general, the 
orbital degeneracy can be lifted by lowering the crystal 
symmetry. For a rough estimation, we assume that the 
orbital degeneracy in pyramidal positions is lifted, which 
gives the following spin-orbital degenaracies: g2o = 12, 
.92p = 4, gzo = 1, gzp = 3, gio = 6, and g^p = 2. 

The hatched area in Fig. 4(a) shows the entire avail- 
able range of A(5„(a:) and A(5p(x) for i?BaCo205+2:, 
which is obtained by changing the equilibrium constant 
K from one extreme to another (where p changes from 
to 1) and keeping the spin-orbital degeneracies un- 
changed. As can be seen in Fig. 4(a), all experimen- 
tal values of the Seebeck coefficient lie within this area. 
This gives confidence that the entropy contribution of 
charge carriers, which includes their spin-orbital degen- 
eracy, can account for the doping dependence of the See- 
beck coefficient and its remarkable divergence at x = 0.5. 
Also, this model naturally accounts for the somewhat en- 
hanced thermoelectric response in electron-doped crys- 
tals in comparison with the hole-doped ones: The larger 
Seebeck coefficient for electrons can be explained by the 
larger degeneracy of HS-Co^+ in comparison with LS- 
Co4+. 

In conclusion, the present study of thermoelectric and 
transport properties in i?BaCo205+2; gives a solid exper- 
imental support to the idea that strong electron correla- 
tions and spin-orbital degeneracy can bring about a large 
thermoelectric power in transition-metal oxides. 
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